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Abstract 

A novel charged particle tracking device, a high precision time projection chamber 
with gas electron multiplier and pad readout, is a leading candidate as the cen- 
tral tracker for an experiment at the International Linear Collider. To characterize 
the performance of such a system, a small TPC has been operated within a mag- 
netic field, measuring cosmic-ray and laser tracks. Good tracking resolution and two 
particle separation, sufficient for a large scale ILC central tracker, are achieved. 
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1 Introduction 



Time projection chambers (TPCs) have been an important part of many large 
particle physics experiments since their initial development in the 1970's. [1,2] 
Traditionally, ionization tracks are imaged at the endplates with wire grids, 
which provide gas amplification, and pads. The signals sensed on anode wires 
and cathode pads are predominantly due to the motion of positive ions away 
from the gas amplification region. When operated in a magnetic field parallel 
to the drift field, the momentum resolution of the device is limited by non-zero 
ExB [3] in the vicinity of the wire grids. The multi-track resolution is limited 
by the wide pad response function and the slow drift velocity of the positive 
ions. 
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A time projection chamber is a leading candidate to be the main tracker for an 
International Linear Collider experiment. [4] In order to improve the momen- 
tum and multi-track resolutions, as required by the physics objectives, it has 
been proposed that the wire grids be replaced by a micropattern gas avalanche 
detector, such as a gas electron multiplier (GEM) [5] or micromegas [6]. For 
these types of devices, ExB can be negligibly small. Furthermore, since the 
pad signals are due to the motion of electrons as they approach and arrive on 
the readout pads, the spatial extent of the signals can be much narrower and 
their risetimes much faster, thereby improving the multi-track resolution. 

In fact, the signals can be so narrow as to present a challenge for conven- 
tional pad readout in a large detector. The outer radius of a linear collider 
experiment main tracker is envisaged to be approximately 2 m and to keep 
electronics costs reasonable, the readout pads need to be no smaller than 
about 2 mm x 6 mm. [4] A strong magnetic field of 4 T is being considered 
to reach the momentum resolution goal. In this magnetic field, gases with fast 
drift velocity at low drift fields can have transverse diffusion contants of about 
30 /im/^/cm. To achieve optimal transverse resolution with 2 mm wide pads, 
a mechanism to defocus the drifting electron charge cloud after amplification 
is required, so that the signals are sampled by at least 2 pads per row. 

One way to defocus the charge cloud is to use gas diffusion between the GEM 
foils and the readout pads. In these regions the electric field can be much 
larger and by selecting an appropriate gas, the diffusion constant can be much 
larger there than in the drift region. This is the method that is investigated 
in the work reported here. Ideally, the defocusing should be sufficient so that 
the standard deviation of the charge clouds should be at least \ — | of the 
pad width, as explained in appendix A. An alternative defocusing method, 
particularly well suited for micromegas devices, is to place a resistive foil 
above the pads. [7] 

The remaining sections of the paper describe studies with a TPC consisting 
of a drift volume coupled to a double GEM structure and conventional pads. 
In section 2, the design and general properties of the TPC and the laser 
delivery system are described. The data samples collected with this device are 
summarized in section 3, and the general properties of the data are shown. 
Section 4 describes a program used to simulate cosmic and laser tracks in 
the TPC, which was used to develop and test reconstruction algorithms. In 
section 5 the methods employed to fit tracks of ionization in the TPC are 
described in detail. Section 6 shows the performance of the TPC, as measured 
under a variety of operating conditions, and compares these to results from 
simulations. 
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Fig. 1. Design drawings of the readout module (left) and drift module (right). The 
two modules are inserted into either ends of a cylindrical acrylic tube to make a 
closed gas volume. The readout module shows the 8 front end cards of the electronics 
readout system inserted. 



2 TPC and laser system 



After gaining experience operating a 15 cm drift length TPC with GEM read- 
out without magnetic fields [8], a new TPC was designed with a 30 cm drift 
length, specifically to be deployed in magnets available for use at the TRIUMF 
and DESY laboratories in 2003. Following successful cosmic ray tests in the 
magnets that year, a laser delivery system was constructed for further tests 
with the TPC in the DESY magnet in 2004. 



2.1 TPC design 



The TPC is a cylindrical design, constructed primarily with acrylic. The drift 
length is 30 cm with a field cage of brass hoops with a pitch of 5 mm held in 
slots in an acrylic holder. The potentials of the hoops are provided by a chain 
of 10 MQ surface mount resistors in the gas volume. The cathode endpiece is a 
circular piece of G10 with copper cladding. The other end of the drift volume 
is terminated by a circular endpiece with a 10 cm square hole, equal to the 
size of the GEM foils. The endpiece has copper cladding, and narrow wires 
strung across the hole, to provide a well defined termination of the electric 
field. The drift volume section of the TPC is inserted into an outer acrylic 
cylinder, with an outer diameter of 8.75 in. A readout module, holding a pad 
array constructed on a printed circuit board and two GEM foils, is inserted 
into the cylinder on the other end, to make a closed gas volume. The readout 
and drift modules design drawings are shown in Fig. 1. 
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Fig. 2. The configuration of the high voltage system is illustrated, (a) The orienta- 
tion of the components, (b) The high voltage for the drift volume is controlled by 
two power supplies (maximum ratings shown). The rectangles show resistances in 
MO. The 100 MO drain resistance is provided so that different drift fields can be 
applied without the lower supply sinking current, (c) The high voltage distribution 
for the GEM foils. The 1001 MO resistances indicate the 1000:1 dividers for moni- 
toring the voltages. Not shown are the isolated voltmeters across shunts to monitor 
the GEM currents. 



2.2 Drift field and GEM operation 



A schematic diagram of the high voltage configuration is shown in Fig. 2. All 
resistors shown are installed in a HV distribution box, except for the field 
cage network, which are surface mount resistors located in the TPC volume. 
The HV distribution box also contained 100 kf2 shunt resistors on one line to 
each GEM, to which isolated digital voltmeters are attached. This allowed for 
current monitoring at the level of 1 nA on each GEM. 

During data taking, electric fields in the drift volume were 90-250 V/cm, 
chosen to maximize the drift velocity, and the field between the drift endpiece 
and GEM 1 was 10 V/cm larger. The transfer field between the two GEMs 
was about 2.5 kV/cm, and the induction field to the readout pads was about 
3.5 kV/cm. The potential across each GEM was 370-380 V, providing an 
effective gas gain of the system in the range 5-8 xlO 3 . 



2.3 Readout electronics 



An array of rectangular readout pads, roughly 10 mm 2 in area are laid out 
on a printed circuit board, arranged in rows of 31 or 32 pads across. The 
pads are at nominal ground potential, and each row is connected to a front 
end card developed for the STAR TPC [9]. The card contains charge sensitive 
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Fig. 3. A schematic drawing (not to scale) of the laser system in plan view. (PD) 
is a photodiode used for triggering the readout electronics. (FS) are fused silica 
windows, used to reflect the beam and reduce its intensity. (BS) are beam splitters. 
(BB) are movable beam blockers to block either beam. The movable mirror allows 
the horizontal separation of the beams to be adjusted. The beams are transported 
below the TPC, and are reflected up through the TPC drift volume between the 
field rings, by a movable mirror, and optionally by a movable beam splitter. 



preamplifiers-shapers and when a trigger is received, the amplified signals are 
sampled and stored into a 512 time bin switched capacitor array. The card 
was modified to increase the pedestal to allow negative going pulses to be 
sampled with a large dynamic range. The signals are digitized asynchronously 
and 10 bits per channel per 50 ns time bin are stored in the raw data files. 
For cosmic ray studies, trigger signals are formed from the coincidence of 
scintillator paddles placed above and below the TPC. For laser studies, trigger 
signals are provided by a photodiode. 

2.4 Laser system 

In order to perform controlled studies of the GEM TPC, a laser delivery 
system was constructed for use in the DESY magnet system with an ultraviolet 
(266 nm wavelength) Nd:YAG laser. The system is able to provide a single 
beam or a pair of beams to the TPC perpendicular to the drift direction, at 
any location along the drift distance, under remote control. The TPC outer 
acrylic tube was fitted with long quartz windows to accept the laser light. A 
schematic drawing is shown in Fig. 3. 

The beam from the laser is reflected by two fused silica slides (each with the 
back face sandblasted to reduce ghost reflections) to reduce the beam intensity 
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and focused by a F = —30 mm plano-concave lens followed by a F = 100 mm 
plano-convex lens. The beam is split, with the reflected beam providing a 
fixed direction, and the through going beam reflected off a movable mirror to 
provide a second beam parallel to the fixed beam. The separation of the two 
beams can be adjusted to between and 12 mm. Beam blockers can move 
into or out of either beam to provide a single beam to the TPC. 

When the TPC is placed in the laser delivery system holder, a movable splitter 
and mirror are underneath the TPC to reflect the beam into the TPC volume 
at two drift distances simultaneously. The splitter can be rotated out of the 
beam so that the light enters at only one drift distance. A steering mirror was 
also in place to allow a transverse separation of beams that enter at two drift 
distances. 

Remote control was provided through the use of servos and stepper motors 
controlled via the serial and parallel ports of the data acquisition computer. 
Operation under remote control was necessary for personal safety considera- 
tions regarding the strong magnetic field and the ultraviolet laser. 



3 Data samples and characteristics 

In 2003, cosmic ray data were collected at a variety of magnetic field strengths 
for two different gas mixtures: P5 (Ar:CH 4 95:5) and a mixture referred to as 
TDR-gas (Ar:CH 4 :C0 2 95:3:2). The pad pitch was 2 mm x 7 mm. Initial 
tests were performed in the 1 T warm magnet at TRIUMF followed by a 
longer data taking period in the superconducting magnet at DESY. Two of 
the electronics cards were faulty during much of the data taking run, but 
6 rows were found to provide sufficient information to characterize the tracking 
resolution performance. 

In 2004, cosmic ray and laser calibration data were collected primarily at a 
magnetic field of 4 T at DESY, with the same two gas mixtures. Two readout 
pad boards were used: the original with 2 mm x 7 mm pads, and a new 
board with 1.2 mm x 7 mm pads. The study with the narrower pads was 
performed because it was found in the 2003 run that neither gas provides 
adequate defocusing between the GEM foils for the 2 mm wide pads at high 
magnetic fields. Gas mixtures with higher concentrations of methane, such 
as P10, could provide the required defocusing at 4 T, but this gas was not 
permitted in the DESY magnet test area. The initial data taking in 2004 was 
with the wider pads, after which the TPC was briefly opened to insert the 
narrow pad plane. The data taking started with P5 gas when it was found 
that the drift velocity at 90 V/cm, where it is expected to be maximum, was 
found to be much lower than expected and increasing with time. By operating 



6 



Table 1 

List of cosmic data samples recorded in 2004, in order of data taking. The name 
assigned to the each configuration is referred to in subsequent tables. The first data 
set likely contained a large concentration of water. 



Name 


Gas 


B field 
[T] 


Pad pitch 
[mm] 


Drift field 

[V/cm] 


p5B4w 


"P5" 


4.0 


2.0 


160 


tdrB4w 


TDR 


4.0 


2.0 


230 


P 5B4n 


P5 


4.0 


1.2 


90 


tdrBln 


TDR 


1.0 


1.2 


230 


tdrBOn 


TDR 


0.0 


1.2 


230 


tdrB4n 


TDR 


4.0 


1.2 


230 



at 160 V/cm, the drift velocity was maximised and found to be stable. After 
the data taking period was completed, it was found that long lengths of the 
gaslines were copper that had been previously exposed to atmosphere, and a 
gas analyzer showed that under these conditions, significant amounts of water 
remain in the system for a few weeks after starting the gas flow. The drift 
velocity and diffusion of TDR gas is less sensitive to water. 

The 2003 data had several faulty electronics channels, unlike the 2004 data, 
and since the laser system was only incorporated for the 2004 data set, the 
focus of the data analysis reported in this paper is on the 2004 data samples, 
summarized in table 1. 

A typical cosmic ray event is illustrated in Fig. 4, along with representative 
pulses on some of the pads. Pads near the center of the track have signals 
due to the collected charge that are significantly different in shape than pads 
further away, which have only signals induced by the motion of electrons in 
the induction gap. Due to the shaping of the readout electronics, the induced 
pulses are transformed into bipolar pulses. 

To estimate the charge collected by pads, clusters of signals in a row are 
searched for, and the time bin with the largest sum of signals in neighbouring 
pads is taken to be the arrival time of the center of the electron cloud. For each 
pad in the row, the sum of the pulse height in that time bin and the three 
previous and three following time bins is formed to estimate the collected 
charge. Studies show that the induced signals roughly cancel out with this 
algorithm. 
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Fig. 4. A typical cosmic ray track observed with the wide pad layout at 4T with 
P5 gas. The brightness of the pads indicates the amount of charged collected. The 
three very wide pads are used during data taking for selecting events that contain a 
track in the active area. The pulses for 4 neighbouring pads are shown on the right. 



Figure 5 shows cosmic ray events at the same drift distance recorded with 
different magnetic fields, and the effect of reduced diffusion in the drift volume 
due to the magnetic field is apparent. 



4 Simulation 

A simulation package was created in order to develop reconstruction methods 
and better understand the various data distributions. [10] Ionization electrons 
can be inserted into a gas volume either along a straight trajectory using 
a parameterization of ionization fluctuations (to simulate muon tracks) or 
Poisson fluctuations (to simulate laser tracks), or along a realistic trajectory 
by using the energy loss of particles traversing the gas, as calculated by the 
GEANT3 [11] program. For most of the simulation studies presented here, the 
latter approach was used with GEANT3 description of the DESY test setup 
and a cosmic muon generator [13], in order to reproduce the cosmic muon 
acceptance. 
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(a) (b) (c) (d) 

Fig. 5. Typical cosmic ray tracks observed with the wide pad layout at different 
magnetic fields with P5 gas. For all events the drift distance is approximately 25 cm. 
The brightness of the pads indicates the amount of charged collected. The black pads 
are disabled due to electronics problems, (a) T, (b) 0.9 T, (c) 2.5 T, (d) 5.3 T. 



The ionization electrons are transported in the gas towards the TPC endplate 
according to a specified drift velocity and diffused according to given transverse 
and longitudinal diffusion constants. The electrons enter the GEM foil holes, 
with a given collection efficiency, are amplified according to either a Poisson 
or an exponential gain function, and exit the holes with a given extraction 
efficiency. Between the foils and between the second foil and the pad plane they 
diffuse, with generally larger diffusion constants. The electrons are collected 
by an array of pads, and the pad signals are shaped, digitized, and saved in 
simulated data files that can be analyzed by the same software that analyzes 
the real data. 

For the simulated samples shown in this paper, the drift velocity and dif- 
fusion constants were set to be the values deduced from the data or from 
Magboltz [12] simulations, the collection and extraction efficiencies set to 1, 
and an exponential gain function assumed for each GEM foil. 



5 Data Analysis 

As part of the simulation package [10], methods were developed to find and 
reconstruct ionization tracks. A simple track finder uses the largest clusters 
from the outermost rows to define straight trajectories as a starting point for 
the track fitter. Track fitting is performed separately in the readout pad plane 
and in the vertical-drift plane. The spatial resolution is most critical in the 
readout pad plane, since it directly determines the momentum resolution of a 
TPC. To achieve the ultimate resolution, a special track fitting algorithm was 
developed, and is described in detail below. No channel by channel calibration 
constants or empirical corrections were applied in the analysis. 
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5. 1 Tracking in the pad readout plane 



To accurately reconstruct the ionization tracks from the measurements of 
charge collected by the readout pads, a maximum likelihood track fitting al- 
gorithm was developed. This takes into account the non-linear charge sharing 
function, which is necessary to achieve good resolution when the pads are 
as wide as the distribution of electrons arriving on the pads. An alternative 
method, using a simple linear centroid determination of the track coordinates 
as it crosses each row of pads, was found to perform poorly under these con- 
ditions. 

A right handed coordinate system is used to define a set of track parameters, 
x , 00) l/ r ! an d cr, that describe ionization tracks projected into the pad 
readout plane. The x — y plane is parallel to the pad plane and perpendicular 
to the drift direction, the +y direction is the vertical direction, and the +z 
direction is opposite to the drift direction. The choice of the coordinate origin, 
(x = 0, y = 0, z = 0), is chosen to be at the centre of the readout pad array. 
The location of the track as it crosses the x axis is x and its local azimuthal 
angle at that location is O - Vertical tracks are assigned an angle O — 0, and 
positive angles correspond to right handed rotations of the track about the z 
direction. To account for the helical trajectories of charged particles in an axial 
magnetic field, the inverse radius of curvature, l/r, is included, with positive 
values assigned to those whose centre of curvature is displaced with respect 
to the track in the +x direction. The final track parameter, a, is the standard 
deviation of the charge distribution about the trajectory as measured at the 
pads. 

The maximum likelihood track fit is described in detail in appendix B. A 
relatively simple model is used to describe the expected charge distribution 
across the readout pads, given the track parameters defined above. As usual, 
estimates for the track parameters for each event are made by maximizing 
the likelihood function. Provided enough rows of data are recorded and good 
starting values for the parameters are used, the maximum of the likelihood 
function is converged upon quickly using standard numerical techniques. 

The second derivatives of the likelihood function can be used to estimate the 
covariance of the track parameters. The error matrix elements estimated in 
this fashion are sensitive to the model assumptions and its parameters, and 
therefore the resolution are estimated directly from the data, by measuring 
the standard deviation of residuals from repeated measurements. 

One effect that causes the errors estimated by the fit to be underestimated is 
the non-uniform ionization that occurs within a row. Except for tracks that 
are parallel to pad boundaries, this increases the variance of the probability 
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density function for charge deposition on pads near the track. The resulting 
degradation of track resolution for non-vertical tracks is referred to as the 
"track angle effect". Another important effect that degrades the resolution is 
the variance of the single-electron response of the GEM, which reduces the 
statistical power of the original ionization electrons. 

One effect that causes the errors to be overestimated is the model assumption 
that the diffusion within and after the GEM structure is small. The observed 
distribution of charge across the pads is used to estimate the diffusion of the 
ionization electrons, whereas in reality, the distribution is due both to the 
diffusion in the drift volume and the diffusion of the amplified electrons in the 
GEM structure. Since the multinomial distribution that is used refers to the 
original ionization electrons, the estimated uncertainties from the track fit will 
be overestimated. 

5.2 Two track fitting in the pad plane 

If two ionization tracks are near enough, some pads will collect charge from 
both tracks, and the tracking resolution for each track will degrade. To analyze 
data with two nearby tracks, the maximum likelihood track fit was modified 
to include the charge from two tracks. To allow for non-uniform ionization 
along the track direction, the ratio of the charges sampled from the two tracks 
is treated as an additional nuisance parameter for each row. In other words, 
the relative amplitudes of the pulses from each track are not assumed to be 
one half, or some other fixed value, but instead the relative amplitudes for 
each row are varied in the maximization of the likelihood. In these studies, 
the width of the charge distribution for each track, a, is fixed. 

5.3 Tracking in the vertical- drift plane 

In the vertical-drift plane, a more traditional approach is applied. For each 
row, the drift time is determined and converted to drift distance. The set 
of points defined by the pad row centres and the drift distances are fit to a 
straight line, by minimizing the squares of the deviations. 

5.4 Selection criteria 

A study of noise present in all channels identified three bad channels for the 
entire 2004 data run and a few single bad channels for subsets of the data run. 
The information from these pads is ignored in the analyses presented here. 
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Table 2 

Measurements of the drift velocity, diffusion, and defocusing, for the 2004 cosmic 
data sets compared to expected values. The data set names are defined in table 1. 
The columns labelled "vj sim" and "D sim" shows the results from a Magboltz [12] 
simulation, whereas the final column shows the results from the full TPC simulation, 
using Magboltz values as the input parameters for the gas properties. 



Data 


Vd 


Vd sim 


D 


D sim 


O"0 


do sim 




[cm///s] 


[cm/H 


[/Ltm/-y/cmJ 


\pm/y/cm\ 


H 


H 


p5B4w 


3.84 ± 0.08 


3.64 


76 ±5 


67 ± 1 


429 ±2 


350 ±2 


p5B4n 


3.85 ± 0.04 


4.14 


34 ±5 


43 ± 1 


382 ± 1 


369 ± 1 


tdrB4w 


4.51 ±0.05 


4.52 


71 ± 10 


69 ± 1 


367 ±4 


262 ± 1 


tdrB4n 


4.54 ± 0.06 


4.52 


70 ± 5 


69 ± 1 


319 ±3 


255 ± 1 


tdrBln 


4.66 ± 0.06 


4.52 


205 ± 10 


206 ±2 


509 ±2 


289 ±2 


tdrBOn 


4.68 ± 0.06 


4.52 


348 ± 20 


468 ± 10 


918 ± 15 


580 ± 1 



Events are selected containing a cluster of signals in each row, with the thresh- 
old set to be typically over 99% efficient per row for cosmic ray tracks. Multiple 
track events are removed by requiring no more than one cluster per row. The 
reconstructed track is required to have its x$ coordinate well within the ge- 
ometric acceptance of the pads, and the azimuthal angle is required to be 
1 0o I < 0.1. Likewise, the z coordinate is required to be well within the drift 
volume, and the dip angle A is required satisfy, | tanA| < 0.6. A small number 
of events are removed by excluding those with an anomalously large a which 
arise when a very large ionization, due to a delta ray, cause large induced sig- 
nals on a large number of pads in a row. Finally a few events are removed by 
requiring that the track fit converge properly, and the estimated errors on xq 
and a be within a reasonable range, to ensure that the fit results are sensible. 



6 Results 

6.1 General properties 

The drift velocity is quickly determined with the laser system by measuring 
the arrival times of laser tracks at known drift distances. This was also used to 
ensure that the drift velocity was relatively constant throughout a data run. 
The results are shown in table 2. 
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Fig. 6. On the left, the points show the variance (ie. a 2 in mm 2 ) and the zq for 
events in TDR gas at B = 1 and 4 T. The right plot shows this data divided into 
10 bins of drift distance, with the mean and error in the mean shown, along with 
linear fits to the data. 



The diffusion properties of the gas, both in the drift volume and in the ampli- 
fication section determine, to a large extent, the performance of the TPC. The 
track parameter, a, gives an indication of the diffusion present for each event. 
Figure 6 shows the variance (a 2 ) for cosmic tracks at various drift distances 
in TDR gas at B = 1 and 4 T. The dependence is linear, as expected from 
diffusion. 

An estimate of the diffusion constant, D, is given by the square root of the 
slope and the defocusing, a , is given by the square root of the intercept. 
When the technique is applied to simulated data, for a broad range of input 
diffusion constants, the estimated diffusion constants are found to be about 
10% too small, when the noise parameter, described in appendix B, p no i se is 
set to 0.01. Furthermore, the discrepancy increases for larger noise parameter 
values. The simulation, however, appears to describe the effect, since the ratio 
of diffusion constants in data to those in the simulated data is found to be 
constant within 3% for 0.001 < p no is C < 0.05. Table 2 shows the results for 
diffusion and defocusing as measured by the cosmic ray data collected in 2004, 
compared to the expectation. The measured diffusion constants have been 
corrected by scaling by a factor of 1.11 ± 0.03. 

The observed drift velocity and diffusion constants are in relatively good agree- 
ment, apart from the diffusion at zero magnetic field. The defocusing term is 
larger in the data as compared to the simulation, for all data sets. 

The energy loss of cosmic muons traversing the gas is measured using all 11 
pad rows, including the three very wide pads. Shown in Fig. 7 is the truncated 
mean number of electrons collected on a pad row per mm of path length 
sampled by the row, where the single largest and smallest values are discarded 
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Fig. 7. The left figure shows the dE/dx distribution in terms of the truncated 
average number of electrons collected on a pad row per mm of path length sampled 
by the row. The points show the data from the p5B4w sample where the GEM 
system gain is estimated to be approximately 8 x 10 3 for this run. The shaded 
histogram shows the result from the simulation. The right figure shows the mean 
dE/dx plotted in bins of momentum for data and simulation. 



in calculating the truncated mean. It is seen that the GEANT3 simulation 
gives a somewhat broader distribution than observed in the data. The increase 
in dE/dx as a function of estimated momentum, as deduced by the track 
curvature, is seen, as expected. Overall, the dE/dx resolution is approximately 
17%. In comparison, an empirical relation [14] that reproduces the dE/dx 
performance of many large scale experiments predicts 16% for this setup. The 
modest resolution is a result of the short overall sampling length of 86 mm. 



6.2 Track resolution in the pad plane 



The likelihood track fit procedure does not involve individual space points 
along the length of the track, so there is no possibility to use the scatter of 
such points to estimate the transverse resolution per point. Instead, the reso- 
lution is deduced by comparing the transverse coordinate estimates from fits 
using different sets of rows. First, the likelihood fit to all 8 rows of data is 
used to define a reference. Next, the data from a single row is used to esti- 
mate the horizontal coordinate, in a fit that leaves the other track parameters 
fixed according to the reference. This is done for all events in a data set and 
the distribution of residuals, between the reference and the single row fit, is 
fit to a single Gaussian to determine the standard deviation, a w . In general, 
the residual distributions are found to be well described by a single Gaus- 
sian. Because the data from the single row was used in the reference, o~ w will 
underestimate the true resolution. To overcome this, the process is repeated, 
but this time the reference fit does not use the information from the single 
row. The corresponding standard deviation, <t wo , is an overestimate of the 
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Fig. 8. The spatial resolution in the pad plane for a row of pads with 7 mm height 
is shown for P5 (left) and TDR (right) gases as a function of drift distance. Wide 
(2 mm) and narrow (1.2 mm) pitch pads are compared, along with results from 
simulated data sets (open points, labelled MC). 



true resolution because the reference track is not perfectly measured itself. 
To properly estimate the resolution from a single row, the geometric mean of 
these two standard deviations is used. This approach is exact for a traditional 
least squares fit [8] , and is also found to work well in simulations and with the 
laser data, as described in Appendix C. 

The single row resolution for cosmic ray tracks measured in a 4 T magnetic field 
by the TPC is shown in Fig. 8 as a function of drift distance and compared to 
the results from simulated data sets. For both gases, the resolution is improved 
substantially by reducing the pad pitch from 2 mm to 1.2 mm, indicating that 
the limited charge sharing for the wide pads is a significant factor in the 
resolution for those pads. Likewise, better resolution is seen in P5 gas than 
TDR gas, which would be expected, given that the defocusing is larger in the 
P5 gas, thereby providing more charge sharing. The overall resolutions, by 
combining events from all drift distances are summarized in table 3. 

The resolution is substantially poorer at lower magnetic fields as shown in 
Fig. 9. A stronger dependence on drift distance is seen as expected, due to the 
increased transverse diffusion constant for the gas at lower fields. 

Figure 10 shows that the resolution degrades with increasing azimuthal angle, 
as expected from the non-uniform ionization of cosmic rays. [8] The depen- 
dence is stronger for the narrow pads than wide pads. The simulation has a 
stronger dependence than seen in the data, which could be related to the fact 
that the ionization fluctuations appear to be larger in the simulation sample, 
as shown in Fig. 7. 

The dependence of the location that a track crosses a pad on the bias and 
resolution for a pad row, is shown in Fig. 11. In the figures, the local pad 
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Table 3 

Overall single row spatial resolution in the pad plane for the 4 data sets collected 
at 4 T and the corresponding simulated data sets. 



dataset 


Resolution 
\pm] (data) 


Resolution 
[/im] (sim.) 


p5B4w 


108 ±1 


92 ± 1 


p5B4n 


68 ±1 


68 ± 1 


tdrB4w 


117 ±2 


100 ±1 


tdrB4n 


83 ±1 


87 ± 1 




10 15 20 

drift distance (cm) 



25 30 



Fig. 9. The spatial resolution in the pad plane for a row of pads with 7 mm height 
is shown for different magnetic field strengths, for the TDR gas with narrow pads. 



coordinate is the location of the reference track (fit from all rows) as it crosses 
the pad, represented as a fraction of the total pad width. A negative bias is 
seen for tracks with a somewhat negative local pad coordinate and a positive 
bias is seen for tracks with a somewhat positive local pad coordinate. This is a 
consequence of the fact that the width of the charge cloud, a, is underestimated 
because of the use of the noise parameter, p no ise, as described in section 6.1. For 
smaller values of p n oise ; c increases and the amplitude of the bias oscillation 
decreases. The magnitude of the bias oscillation is larger in the simulated 
samples. In all cases, the bias is smallest at the edges or the centre of the pad, as 
expected, since a does not influence the reconstructed transverse coordinate at 
those locations. Overall, the effect is small enough that it does not significantly 
contribute to the overall pad row resolution. 
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Fig. 10. The spatial resolution in the pad plane for a row of pads with 7 mm height 
is shown for different azimuthal angles is shown for P5 data samples at 4 T. 
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Fig. 11. The single row bias and resolution for P5 data at 4T is shown as a function 
of local pad coordinate (the location of the reference track on the pad, mapped to 
—0.5 - 0.5, where the centre of the pad is located at 0). The left figure shows a bias 
dependence that arises from an underestimate of the charge width, a. The right 
figure shows that for the wide pads, the resolution is better near the pad edges. 



For the wide pads, the resolution is best for tracks that cross near the edge 
of the pad. This dependence is expected and described in Appendix A. The 
simulation shows the same trend, but shows a much more peaked behaviour at 
the centre of the pad. For the narrow pads, the resolution is roughly constant 
across the pad for data and in the simulation. 
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narrow, z = 54 mm 
narrow, z = 263 mm 
laser simulation 
laser simulation 
laser simulation 
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3 4 
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Fig. 12. The degradation in transverse resolution is shown as a function of the 
transverse separation between two parallel tracks, for P5 gas at 4T. The z values 
shown are the drift distances of the tracks. The laser data and simulations (open 
symbols with the same shape) show rough agreement. The simulation result for 
muons with 2 mm wide pads, indicate good resolving power for separations above 
around 3 mm. 



6.3 Two track resolution in the pad plane 



To study the capability of the detector to extract track coordinate information 
for nearby tracks, two parallel laser beams were brought close together at the 
same drift distance. This represents the most difficult situation for resolving 
two tracks in a TPC. By using the beam blockers, events were recorded with 
individual beams, as well as events with both beams present. As described in 
section 5.2, the likelihood track fitter was modified to fit two tracks for this 
study. The fractional increase in the standard deviation of the Xq coordinates 
when both beams are present is used to quantify the two track separation 
power. A simulation of the laser tracks was also performed, where the laser 
ionization was assumed to be Poisson distributed. Figure 12 shows the degra- 
dation in the resolution as the beams are brought together for the data and 
simulation, in wide and narrow pads in P5 gas at 4T. There is rough agreement 
with the data and simulation. To estimate the two track resolution capability 
for minimum ionizing tracks in the device, a simulated data set of muon pairs 
at different separations was generated for the wide pads, and the results are 
included in Fig. 12. The track fitter fails to converge for 1% (5%) of the muon 
pair events for track separations of 3 mm (2 mm). In summary, this study 
indicates that for a GEM TPC with 2 mm wide pads, significant information 
can be extracted from a pad row when two tracks are as near as 3 mm in that 
row. 
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Table 4 

Overall single row spatial resolution in the drift direction for the 2004 data sets 
and the corresponding simulated data sets. 



KA. IA> Lj CIO V_ V 


Resolution 

[mm] (data) 


Resolution 
[mm] (sim.) 


p5B4w 


0.99 ±0.01 


0.83 ± 0.01 


p5B4n 


1.02 ±0.01 


0.81 ±0.01 


tdrB4w 


0.99 ±0.01 


0.85 ±0.01 


tdrB4n 


0.95 ±0.01 


0.87 ±0.01 


tdrBln 


0.99 ±0.01 


0.90 ± 0.01 


tdrBOn 


1.15 ±0.02 


0.93 ± 0.01 



6.4 Track resolution in the drift direction 

The focus of this study is on the resolution performance in the pad plane for 
which the demands are the greatest and the magnetic field has a direct impact. 
For completeness, the observed resolutions in the drift direction are reported in 
this section. A linear track fit is performed using the vertical coordinates of the 
pad rows and the drift time for the cluster observed on the row. The residuals 
of the drift time from each pad row with respect to the full track fit are fit 
to Gaussians to determine the standard deviations. As before, the geometric 
mean of the standard deviations determined when the track fit includes or 
excludes the pad row is taken to be the resolution. The resolution is seen to 
degrade with drift distance, as expected from the longitudinal diffusion. The 
data resolution is somewhat poorer than the values found from the simulation, 
as shown in table 4. Although the TDR mixture has a lower longitudinal 
diffusion constant than P5, it has a faster drift velocity, and with the 50 ns 
sampling, the resolutions of the two gases are found to be quite similar. 



7 Conclusions 

The results presented in this paper indicate that a GEM TPC with pad read- 
out is a viable candidate for the central tracker for an experiment at the 
International Linear Collider. When operated in a 4 T magnetic field with 
modest size pads, a full size TPC could collect some 200 pad-row measure- 
ments each with transverse spatial resolution of approximately 100 /im, which 
should be sufficient to achieve the transverse momentum resolution goal of 
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a(l/p t ) = 2 x 10~ 4 (GeV/c)" 1 . [4] In addition, the resolution does not de- 
grade significantly for nearby tracks, provided their separation is more than 
about 1.5 times the pad width. 

The relatively simple simulation package reproduces many of the general fea- 
tures seen in the data. It should be useful, therefore, in optimizing the config- 
uration of GEM TPCs for future experiments at the ILC and elsewhere. 
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A Resolution and Charge Sharing in a GEM TPC 

Consider a charge cloud of n p electrons arriving at the GEM plane. Due to 
diffusion in the drift volume, the electrons are distributed in the transverse 
direction, x, with standard deviation o^. If the GEMs provide a gain g (with 
negligible variance), the defocusing adds a variance <Tq to the x distribution, 
and if the x coordinate of all of the resulting electrons are measured with no 
uncertainty, then the variance of the mean x coordinate is 



To achieve the diffusion limit, the GEM term (the second term) must be 
much smaller than the diffusion term. In reality, the GEM term has additional 
contributions from the variance in the gain, and the uncertainties in measuring 
the x coordinates for each electron. 




(A.l) 
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In a GEM-TPC the x coordinate for each electron is not measured, but rather 
the charge shared by neighbouring pads is used to deduce the mean x coordi- 
nate. The degradation in resolution due to using large pads can be understood 
analytically. Consider two neighbouring semi-infinite pads with the boundary 
at x — 0. If the electrons are distributed according to the pdf G(x), the ex- 
pectation for the fraction of electrons over the positive-x pad is 

oo 

(F) = J G{x) dx . (A.2) 
o 



If G(x) is Gaussian, with mean /x, standard deviation a, the estimate ft deter- 
mined from the observed fraction F has variance, 

a 2 = 2na 2 e^^a 2 F (A.3) 



where the variance of F is binomial, a 2 F = (F) (1 — (F)) /n and n is the 
number of electrons. If fi = 0, so that the mean of the pdf is at the border 
between the pads, the variance on the estimate ft, is <r| ~ 1.6a 2 /n or in other 
words a factor 1.6 larger than the variance that would result in perfect x 
coordinate measurements for each electron. As the mean of the pdf moves 
away from the pad boundary, the variance gets larger; for fi = la (2a) the 
factor increases to 2.3 (9.0). To keep the GEM contribution to the resolution 
small, the pad width needs to be less than about 3-4 a, so that /x is less 
than 1.5-2 a. The effects of noise and thresholds may limit the pad sizes even 
further. 



B Maximum Likelihood Track Fit in the Pad Plane 



The maximum likelihood track fit uses a simple model to describe the way 
electrons are spread about the readout pads. Due to diffusion, the electrons 
are assumed to be distributed in a guassian fashion about the projection of the 
particle trajectory onto the readout pad plane. The resulting "Line-Gaussian" 
density function is a convolution of the "trajectory density function" and a 
two dimensional isotropic Gaussian probability density function. The standard 
deviation of the Gaussian, a, is a free parameter in the fit, like the other track 
parameters. The fitted value provides an estimate of the diffusion for each 
event. For simplicity in the analysis, a is assumed to be constant along the 
the length of the track, a good approximation for the data samples described 
in this paper. 

The implementation of this model assumes that rectangular pads of equal 
height are arranged in a parallel fashion into horizontal rows and that tracks 
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cross several rows. The height of the row is assumed to be much smaller than 
the radius of curvature, so that the tracks can be considered to be straight 
within the vertical extent of the row. 

The true trajectory density function (tdf) for a typical track is non-uniform 
along its path due to ionization fluctuations. The model, however, assumes 
that within a row, the tdf is uniform. The likelihood, C row , of observing the 
distribution of charge is calculated for each row, and the product of the like- 
lihoods for all the rows defines the overall likelihood function, C In this way, 
the model implicitly allows the tdf to vary from row to row. 

The expected charge collected by a rectangular pad is determined by integrat- 
ing a uniform Line-Gaussian density function over the physical region of the 
pad and is proportional to: 

w/2 h/2 

f f dxdy ( [(x-6)cos0 + ?/sin0] 2 \ 
I(b,<f>,a,h,w)= j j ^exp(^ — j 



rj(b, 0, a, h, w) - rj(b, 0, a, -h, w) 
+rf(b, 0, a, —h, -w) - rj(b, 0, a, h, -w) 



(B.l) 



r}(b, 0,tT, h,w) 



1 ■ ^ \( b+ ^) cos0+ ^sin0,a 

cos (b sm d> \ \ 1 J I , 



. . . u , I u \ a / —u' 

«"' <7) = 2 erf l7sj + 7S exp U^ 



where b is the horizontal distance between the pad centre and the track, is 
the local azimuthal angle of the straight line segment in the pad row, and h 
and w are the height and width of the pads. The local azimuthal angle for a 
row centred at y — y row is 

sin = sin O - y TO w/ r . (B.2) 



The horizontal distance to the track from the pad centred at x — x pa( j in that 
row, can be expressed in a power series of 1/r as 

b = x - x pad - y IOW tan O + \y 2 mw sec 3 O ^ - |y r 3 ow tan O sec 4 O ^ .(B.3) 
In order to calculate the likelihood of the observed charge distribution in a 
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single row, an approximate model is used. The model is developed on the basis 
that signals in a row result from a relatively small number, n pa 0(100), of 
electrons liberated in the ionization process that are subsequently amplified 
in the GEM structure. As an approximation, the distribution of charge is 
assumed to be multinomial with n observations. This approximation is exact 
in the limit that the diffusion within and after the GEM structure is small 
compared to the diffusion in the drift region, the fluctuations due to electronics 
noise is small compared to that due to primary electron statistics, and that 
thresholds are small. This model was initially developed to study data from 
tests without magnetic fields, in which the dominant contribution to diffusion 
was in the drift volume. For certain gases and operating conditions in strong 
magnetic fields, the situation is reversed, whereby the diffusion primarily arises 
in the GEM structure. In fact this feature is desirable in order to ensure 
that there is good charge sharing across relatively large pads, as described in 
section 1. Nevertheless, the model performs very well in either circumstance. 

If Ni electrons are collected by pad % for a system with total GEM gain, g, the 
equivalent number of the original ionization electrons associated to the pad is 
simply defined as 

n, = Ni/g . (B.4) 



These values are non-integer, but a direct continuum extension of the multi- 
nomial distribution, yields the log likelihood function, 

log £ row = ^2 n j log pj + constant (B.5) 

i 



where Pi is the probability for a primary electron from the track to be associ- 
ated to pad i, 

Pi = I (pi, fa, a, h u Wi)/ ^ <j)j, a, h j: Wj) (B.6) 

j 



where j runs over all pads in the row containing pad i. 

A spurious signal in a pad far from the track can cause problems for the like- 
lihood calculation, since the calculated probability for electrons to be present 
there can be vanishingly small. To make the track fit robust to spurious sig- 
nals, the probability for a primary electron to be associated with a pad is 
modified by adding a small constant, p no i S e, 

Pi + Pnoisc /-p, _\ 

Pi -»• T— (B.7) 

1 ~r ^pad Pnoisc 
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where n pa d is the number of pads in the row. The fitted track parameter values, 
apart from a, are found not to strongly depend on the choice of p no ise> f° r the 
data analyzed in this paper, where p no i SC ~ 0.01. 



C Study of resolution determination method 

The method used in this paper to form the estimated single row resolution 
(ESSR) has been tested using simulated samples and laser data. 

For simulated samples, the true values of the track parameters are known. 
Events were simulated with the ionization following a straight trajectory, and 
the resulting pad signals were fit to straight tracks (fixing 1/r = 0). The 
resolution is directly determined from the standard deviation of the residuals 
of the fit xq values with respect to the true values. The direct single row 
resolution is defined to be times this standard deviation, since 8 rows are 
used in the fit. The ESRR is found to agree with the direct result to within 
10%. When the same data is fit to curved tracks (allowing 1/r to float), the 
ESRR is unaffected, but the standard deviation of the residuals for the direct 
result increases due to the correlation between x Q and 1/r. 

Laser data was used in a similar way to check the ESRR. It was found that 
the laser system was quite steady; the mean of the x from track fits to the 
laser data typically drifted by less than 10 fim over a period of 12 hours. In 
the data with two laser beams present at large drift separations, however, the 
reconstructed xq values for the two tracks were found to be correlated. This 
is interpretted as due to an angular jitter of the laser and a simple model is 
used to describe this situation. The measured track parameter x for track i 
(i — 1 or 2 - the near or far track), is described as an outcome of a random 
variable X iy 

Xi = m + R e + Ri (C.l) 

where Ri is a random variable of mean zero describing the jitter of the laser 
system and Ri is a random variable of mean zero describing the random nature 
of the Xq measurement by the TPC. The covariance, ({Xi — ^i)(X 2 — ^2)), 
isolates the contribution from the laser jitter, and its square root is found to 
be about 10 /im, corresponding to an angular jitter of about 5 /xrad. 

Only data with two laser pulses are used, so that the laser jitter term can 
be removed, and the standard deviation of the remainder is multiplied by 
to form the direct single row resolution. The direct single row resolution is 
determined from fits with 1/r = 0, in order that the correlation between xq 
and 1/r not influence the resolution measurement. 
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0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 
direct resolution (mm) 

Fig. C.l. The estimated single row resolution is compared to a direct resolution 
estimate from laser data collected with narrow pads and P5 gas at 4T. The error bars 
indicate the statistical uncertainties of the resolution estimates. In each case, the 
laser jitter contribution has been removed in the calculation of the direct resolution. 
The lines correspond to ESRR = c x direct. 



This is compared to the ESRR for the same data sets in Fig. C.l. The fits 
used in the ESRR calculation allowed 1/r to float, as is done when estimating 
the resolution with the cosmic ray data. The ESRR is found to be somewhat 
smaller than the direct resolution; it appears that ESRR ~ (0.9±0.1)x direct. 

In summary, the simulation and laser studies show that the method used in 
this paper to estimate the single row resolution of a GEM TPC appears to be 
correct within a systematic uncertainty of about 10%. 
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